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Abstract. Advantages introduced by Raman lidar systems for 
cloud base determination during precipitating periods are explored 
using two case studies of light rain and virga conditions. A com- 
bination of the Raman lidar derived profiles of water vapor mixing 
ratio and aerosol scattering ratio, together with the Raman scat- 
tered signals from liquid drops, can minimize or even eliminate 
some of the problems associated with cloud boundary detection 
using elastic backscatter lidars. 


1. Introduction 

Accurate characterization of cloud boundaries is fundamental to 
understanding the radiative effect of clouds on climate. Cloud- 
radiation interactions impact the radiation budget of the Earth, 
both at the surface and in the atmosphere, and are the largest 
source of radiation budget uncertainty [IPCC, 1995]. These inter- 
actions depend on many variables, including cloud depth, cloud 
top height, cloud liquid water content, water vapor profile, and 
solar zenith angle. The calculation of cloud optical depth, for ex- 
ample, requires accurate determination of cloud base and cloud 
top heights. Recent studies suggest that variations in cloud 
boundaries are even more important than variations in cloud mi- 
crostructure for the calculation of cloud optical depth [Kogan and 
Kogan, 1998]. Han and Ellingson [1997] also indicate that cloud 
geometry is the primary factor that controls effective cloud frac- 
tion and is a major uncertainty in calculation of downwelling 
longwave flux at the surface. 

Measurements of cloud horizontal extent predominantly make 
use of satellites [Rossow and Gardner, 1993], while characteriza- 
tion of cloud vertical distributions primarily rely on radar and lidar 
observations [Clothiaux et al. 1998, Uttal et al. 1995]. However, a 
number of difficulties exist in the radar-lidar observation of cloud 
base height. The detected cloud base height is strongly influenced 
by the sensing wavelength as well as the method of detection. 
Occurrence of rain, virga, and radar-confusing bugs in the bound- 
ary layer also complicate cloud base detection. 

This work discusses advantages introduced by Raman lidar in 
the determination of cloud base height using the NASA Goddard 
Space Flight Center's Scanning Raman Lidar (SRL). The SRL 
provides simultaneous profile measurements of aerosol scattering 
ratio (ASR), water vapor mixing ratio (W in g kg 1 ), and relative 
humidity (RH in percent derived using co-located temperature 
profiles from sondes or other sensors). The Raman lidar meas- 
urements will be discussed in relation to the Micro Pulse Lidar 
(MPL), Belfort Laser Ceilometer (BLC) and Millimeter- 
Wavelength Cloud Radar (MMCR). Details of the SRL have been 
reported by Whiteman and Melfi (1999). 
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2. Cloud Base Detection 

Cloud-base detectors, available commercially and widely used 
by the U.S. National Weather Service, like the rotating beam cei- 
lometer and the laser ceilometer, assume cloud base height to be 
the location of the peak in the lidar backscatter signal [Eberhard, 
1986]. This cloud base height definition assumes that the signal is 
small below cloud, rises to a maximum at the start of the cloud, 
and decreases thereafter due to attenuation. Thus, most tech- 
niques used for cloud base height detection search for zero cross- 
ings in the derivative of the returned signal [Pal et al. 1992, 
Campbell et al. 1998, Gaumet et al. 1998]. Other methods look 
for a ’’threshold" value of signal shift from the background [Platt 
et al. 1994] or construct lidar clear-sky power return profiles from 
archived data to test for the presence of clouds [Clothiaux et al. 
1998]. Noise estimates and/or smoothing are used in some of 
these techniques to reject false peaks resulting from cloud inho- 
mogeneity, aerosols and signal noise. 

The above definitions of cloud base height are not as straight- 
forward as they might appear. They are influenced by the instru- 
ment wavelength, and the fact that appropriate definition of cloud 
base height may depend upon the particular research problem to 
be addressed. For surface radiation budget studies, for example, 
cloud base height may refer to the radiatively important part of the 
cloud and neglect any virga-associated excursions, while for visi- 
bility applications, such as pilot/flight safety, alternative criteria 
may be required. 

Although cloud base height is routinely assumed to be the lo- 
cation of the peak in the lidar backscatter signal, it is well known 
that this is an overestimation. The difference in altitude between 
the point where the signal starts to increase and its peak can be as 
large as 800 m [Pal et al. 1992]. Pilot-reported cloud base height, 
for example, is usually lower by as much as 240 m than that de- 
tected by the rotating beam ceilometer [Pal et al. 1992]. The dif- 
ference between the two levels, i.e., where the signal starts to in- 
crease and where it peaks, is a function of the physical interaction 
between the cloud and the instrument signal. Instruments operat- 
ing at shorter wavelengths (near visible) respond to smaller size 
droplets, like those found near the base of non-precipitating 
clouds, that might not be detected at longer wavelengths such as 
used by radar. 

Additional complications arise in the detection of cloud base 
height for high altitude clouds. High cirrus clouds often exhibit 
complex inhomogenous structures and may have rapid changes in 
cloud base height and the lidar signals often have low signal-to- 
noise-ratio as compared to boundary layer clouds. Thus, cloud 
base height algorithms must distinguish between signal shifts due 
to random noise and cloud base height. These differences can be 
challenging in the case of thin cirrus clouds [Sassen and Cho, 
1992]. 

Problems in cloud base height detection also arise when the li- 
dar signal does not undergo a sudden increase (as in fog or the 
initial phase of cloud formation), when multiple peaks occur in 
close proximity (multiple cloud layers not well separated) and at 
times of precipitation or haze between cloud and the instrument. 
Rain and virga introduce large errors in the detection of cloud 
base height, although, in principle, one can differentiate between 
regions of precipitation and cloud from the rate of the slope of the 
backscattered signal [Pal et al. 1992]. 



3. Raman Lidar Contribution 


In addition to the elastic backscatter signal at the transmitted 
wavelength, Raman lidar systems detect signals from water vapor, 
nitrogen and oxygen molecules and from cloud liquid water, due 
to Raman scattering at wavelengths shifted from the laser wave- 
length [Melfi, 1972, Melfi et al. 1997]. Raman lidar measurement 
of water vapor mixing ratio is made using the ratio of signals at 
the shifted water vapor and nitrogen frequencies. The extra in- 
formation provided by Raman lidars, and particularly the water 
vapor and cloud liquid water signals, can be used to resolve some 
of the problems in cloud base height detection. In particular, rela- 
tive humidity calculated from the water vapor mixing ratio, the 
Raman signal from liquid cloud drops, and the attenuation of the 
signals, together with the elastic backscatter at the transmitted 
wavelength, can be used collectively to infer cloud base height. In 
principle, sensitivity to the molecular processes makes the Raman 
signal immune to virga and/or light precipitation and removes the 
dependence of cloud boundary definition on instrument wave- 
length. Thus, information on the water vapor and nitrogen profile 
can be retrieved even during virga and/or precipitation which is 
either not possible or not reliable using other single channel elastic 
lidars. The Raman system is immune to the so-called plankton 
problem, where bugs and other large suspended particles compli- 
cate cloud/non-cloud boundary, and seriously limit radar based 
cloud base height detection accuracy for low level clouds in some 
seasons. 

3.1. 28 September 1997: Virga 

Figure 1 shows analyses of the SRL, Belfort Laser Ceilometer, 
and Micro Pulse Lidar data for 28 September 1997 observed at the 
Department Of Energy Atmospheric Radiation Measurement site 
near Lamont, Oklahoma. Microwave radiometer (MWR) meas- 
ured vertically integrated cloud liquid water path (in mm) are also 
included. Note that the Raman lidar data used here are nighttime 
only. Detailed analysis of wind profiler data, synoptic maps, Ra- 
man lidar and other data sets indicated that the clouds on this day 
were formed as a result of convergent lifting of low level moist 
air. Cloud level winds were less than 10 m sec' 1 . 

As illustrated in Figure 1, cloud base height derived from all 
three instruments (i.e., SRL, MPL and BLC) compared well, ex- 
cept at about 0600, 1300 and 1700. This is not surprising because 
these cloud base heights for all three were derived from the elastic 
backscatter signal (ASR), albeit using different algorithms and 
wavelengths (^j^ vapor ==408nm, ^srl _nitrog« =: 387nm, 

^sRLjaser = 355nm; X BlJ c=9 lOnm; ^Mp L =523nm). The Raman lidar 
cloud base heights were derived using a simple “threshold” 
method where cloud base height is defined as the level at which 
the signal is twice the background value and continues to increase 
(or does not decrease) for at least some distance above. The micro 
pulse lidar cloud base heights were derived using the Scott- 
Spinhime threshold algorithm (Clothiaux et al. 1998). The algo- 
rithm employed by the Belfort Laser Ceilometer, although pro- 
prietary, is similar to that of Pal et al. (1992). 

Around 0600, 1300 and 1700, the lidars generate several differ- 
ent cloud base heights, within a very short time interval, indicating 
problems as they attempt to identify the base height of the cloud 
layers. Coincident observations by MMCR indicated virga condi- 
tions at these times that can confuse elastic backscatter-based 



cloud base height detection. Profiles of SRL (Micro pulse lidar) 
measured aerosol scattering ratio (estimated backscatter), sonde 
and Raman lidar measured water vapor mixing ratio (g kg* 1 ) and 
relative humidity (%) profiles, and the dew point depression (T- 
T d , in °C) at 0531, 0600 and 0830 are plotted in Fig. 2. 

The 0531 sounding (T-T d , dotted line, in Fig. 2A) indicated a 
cloud deck from 2.6 to 3.1 km and a second nearly saturated 
shallow layer from 4.2 to 4.3 km. The lower cloud layer is re- 
solved by Raman lidar-measured aerosol scattering ratio profile 
(heavy dash-dot line). The peak return occurs about 100 m above 
cloud base (twice background). Cloud top height is slightly over- 
estimated by the SRL in comparison to the sonde data. The 
agreement is fairly good. This is also true for the 0830 (Fig. 2D) 
profiles where only the lower cloud layer remains. The 0600 pro- 
files (Fig. 2B), however, are problematic. Although the micro 
pulse and Raman lidar profiles do not correspond exactly (mainly 
due to a difference in signal averaging time and instrument sepa- 
ration), they indicate three or more cloud layers. However, a close 
scrutiny of the data will show that only two layers exist. It is in- 
teresting that a minimum (“dip”) is found about 150 m below the 
0°C level (above a peak at 3.7 km) in both the micro pulse and 
Raman lidar profiles (Fig. 2B). This feature is typical of the so- 
called “lidar dark band” [Sassen and Chen 1995], which is associ- 
ated with the structural collapse and accumulation of melted 
snowflakes and their transition to spherical drops. Thus, the peak 
at 3 km should not be classified as a cloud layer. First, the fact 
that the lidar signals did not get heavily attenuated and gave in- 
formation above the clouds (above 4.2 km) is a good indication 
that the clouds were not optically thick and therefore not deep 
(optical thickness at 0600 between 2.5 to 4.3 km was 1.4, less than 
the 2.6 calculated at 0800 for the layer between 2.5 and 3.5 km). 
Second, the sounding at 0531 (and 0800) indicated sub-saturated 
conditions from 3 km to near 4 km. Finally, the water vapor 
mixing ratio and correspondingly relative humidity values be- 
tween 3 and 4 km far exceed saturation values. This is unrealistic. 
Note the apparent moistening in the layer immediately below the 
lidar dark band that is indicated by the SRL profile at 0600 (Fig. 
2C). This is attributed to additional signal in the water vapor 
channel due to Raman scattering from liquid drops (see Whiteman 
and Melfi, 1999) in addition to evaporation of the falling drops. 
Thus, around 0600, only two cloud layers existed; a layer between 
2.6 and 3 km and another precipitating cloud between 4 to 4.2 km, 
most likely a glaciated altostratus cloud with virga between the 
cloud layers. This is supported by MMCR observations of virga 
around 0600 originating from the upper cloud. In such conditions, 
single channel lidars will have great difficulty in locating the 
cloud base height, even when using sophisticated algorithms. We 
suspect the same complications caused by virga in detecting cloud 
base height by the lidars at 0600 are also the causes for the multi- 
layer cloud indications at 1300 and 1700. 

Excess vapor (a bump in the SRL profile compared to the 
sonde values) was observed near cloud base (2.8 to 3 km) in the 
0830 profile while virga was not reported by the MMCR. Raman 
measurements of water vapor mixing ratio obtained in the pres- 
ence of liquid cloud drops are enhanced by the Raman signal from 
the drops. Separation of the vapor and liquid Raman scattering 
into separate signals will simplify the processes of cloud/virga 
detection. Fortunately, the vapor and liquid water Raman bands 
have sufficient wavelength difference to allow separate detection 
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of liquid and this capability is currently being designed for the 
SRL system [Whiteman and Melfi, 1999]. In the absence of virga 
conditions, the “excess vapor” due to Raman scattering from 
cloud drops is a clear indicator of the onset of cloud water drops 
and hence is a direct indicator of cloud base height. The relative 
humidity, derived from the corrected water vapor profile, will also 
serve as an indirect indicator of cloud base height. In addition to 
cloud base height, the water vapor profiles show the moistening of 
the layer leading to cloud formation, even in conditions of virga 
and precipitation. 

3.2. 27 August 1998: Rain 

Recent improvements in the SRL housing trailer allows for all- 
weather operation, including rain conditions. An example of data 
during such conditions, observed on 27 August 1998 at Andros 
Island, Bahamas, as part of the Convection And Moisture Ex- 
periment (C AMEX-3) field campaign, is shown in Figure 3. 
Aerosol scattering ratio and water vapor mixing ratio profiles are 
shown at two selected times: at .0400 when precipitation reached 
the surface (Fig. 3A) and at 0745 (Fig. 3B). Rawinsonde meas- 
ured dewpoint depression at 0745 is also plotted in Fig. 3C. At 
0400, near-surface aerosol scattering ratio values were about five 
times the background, when precipitation was recorded at the 
ground. Aerosol scattering ratio increased monotonically to 4.2 
km followed by a “dip” to a minimum at 4.6 km and then a peak 
at 5 km, before returning to near-background levels. This peak- 
dip-peak structure near the freezing level is the lidar-dark band 
described previously. A similar shape in aerosol scattering ratio 
and water vapor mixing ratio profiles was also observed at 0745. 
Given only backscatter information (which is the case for most 
elastic lidars), the threshold technique of determining cloud base 
height will not work and most other algorithms would have diffi- 
culties. An alternative solution could utilize available radiosonde 
information or depolarization ratio [Platt et al. 1994]. With Ra- 
man lidar systems, however, the water vapor mixing ratio channel 
may be used. Note that the water vapor channel at 0400 (0745) 
provided usable data up to 5.1 (5.4) km, close to the peak above 
the lidar dark band in the aerosol scattering ratio profile. Above 
5.1 (5.4) km, the lidar signals were severely attenuated (this oc- 
curs for optical depth greater than about 5). The lidar-derived 
relative humidity (using sonde temperature, not shown) indicated 
saturation was reached at an altitude of 5.2 km, slightly lower than 
the peak in aerosol scattering ratio. From the 0745 sonde, the al- 
titude of the smallest dewpoint depression was 5.2 km, in agree- 
ment with the lidar-derived cloud base height, confirming the util- 
ity of the Raman water vapor signal for cloud base height detec- 
tion in rain (0400) or virga (0745) conditions. 


4. Summary 

Problems in current methods of cloud base height determination 
during precipitation conditions using single channel elastic back- 
scatter profiles are discussed. In some cases, the errors can be 
substantial. Principles of a possible new approach, based on ap- 
plication of the NASA/GSFC SRL derived aerosol scattering ra- 
tio, water vapor mixing ratio (and derived humidity profile) and 
cloud liquid water drops (measured as “apparent” excess in water 
vapor mixing ratio) are outlined. Measurement of cloud base 
height using this new approach is demonstrated and applied to 



light rain and non-rain cases. Results of this approach revealed 
that the lidar-based cloud base heights in non-rainy conditions 
compared well with those from micropulse lidar and Belfort laser 
ceilometer. Raman lidar-based retrieval of cloud base height and 
water vapor mixing ratios, during light rain and virga conditions, 
was demonstrated, which is impossible or highly problematic us- 
ing elastic lidars. These results suggest that, where available, 
Raman lidar systems can add substantial information to cloud 
boundary studies in diverse weather conditions. The capability to 
simultaneously utilize the aerosol scattering ratio, water vapor 
mixing ratio, and Raman scattering from liquid water makes Ra- 
man systems a powerful tool for cloud boundary studies. 
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Figure 1. Time series of (a) Microwave Radiometer (MWR) in- 
tegrated liquid water content and (b) cloud base height determined 
from measurements of scanning Raman lidar (SRL, plus), Belfort 
Laser Ceilometer (BLC, open circles), and Micro Pulse Lidar 
(MPL, crosses). Note that in order to make the symbols discerni- 
ble the vertical axis between 2.5 to 3.5 km is expanded and only 1 
in 10 of the points are plotted in (b). The location for the 0 C level 
and the MMCR virga observation times are also indicated. 

Figure 2. Profiles of SRL-measured aerosol scattering ratio 
(ASR, dash-dot) and water vapor mixing ratio (W, solid with error 
bars), and sonde derived dewpoint depression (T-Td, dotted) and 
water vapor mixing ratio (W, thin solid) are plotted at (a) 0531, (b 
and c) 0600 and (d) 0830 UTC. The freezing level is also indi- 
cated. 


Figure 3. Typical profiles of SRL-measured aerosol scattering 
ratio (ASR, dash-dot) and water vapor mixing ratio (W, heavy 
solid) from (a) rain at 0400 UTC and (b) virga at 0745 UTC con- 
ditions for the 27 August, 1998 case. (C) Sonde derived profile of 
the dewpoint depression (T-Td) a * 0745 is plotted for comparison. 
The freezing level, cloud base height and lidar dark-band positions 
are also indicated. 
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